Kinetic models for wealth exchange on directed networks 
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We propose some kinetic models of wealth exchange and investigate their behavior on directed 
networks though numerical simulations. We observe that network topology and directedness yields 
a variety of interesting features in these models. The nature of asset distribution in such directed 
networks show varied results, the degree of asset inequality increased with the degree of disorder in 
the graphs. 
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I. INTRODUCTION 

The distribution of wealth among individuals in an 
economy has been an important area of research in eco- 
nomics, for more than a hundred years [1 0, H 0]. The 
same is true for income distribution in any society. De- 
tailed analysis of the income distribution [1 H[ so far 
indicate 

, , J m a exp(— ra/T) for m < m c , . , 

F{m) ~ \ m-( 1+ ^) for to > m c , (1) 

where P denotes the number density of people with in- 
come or wealth to and a, v denote exponents and T de- 
notes a scaling factor. The power law in income and 
wealth distribution (for to > m c ) is named after Pareto 
and the exponent v is called the Pareto exponent. A 
historical account of Pareto's data and that from recent 
sources can be found in Ref. [1 . The crossover point (m c ) 
is extracted from the crossover from a Gamma distribu- 
tion form to the power law tail. One often fits the region 
below m c to a log-normal form logP(?7i) oc —(log to) 2 . 
Although this form is often preferred by economists, we 
think that the other Gamma distribution form Eqn. ([1]) 
fits better with the data, because of the remarkable fit 
with the Gibbs distribution in Ref. @, 0, H . 

Considerable investigations revealed that the tail of the 
income distribution indeed follows the above mentioned 
behavior and the value of the Pareto exponent v isgener- 
ally seen to vary between 1 and 3 @, E3, El E3, El El Hi- 
lt is also known that typically less than 10% of the pop- 
ulation in any country possesses about 40% of the to- 
tal wealth of that country and they follow the above 
law, while the rest of the low income population, fol- 
low a different distribution which is debated to be either 
Gibbs [3, E3, E3, ES El El or log-normal [H El ■ 

The striking regularities observed in the income distri- 
bution for different countries, have led to several new at- 
tempts at explaining them on theoretical grounds. Much 
of it is from physicists' modeling of economic behavior 
in analogy with large systems of interacting particles, as 



treated, e.g., in the kinetic theory of gases. According to 
physicists, the regular patterns observed in the income 
(and wealth) distribution may be indicative of a natural 
law for the statistical properties of a many-body dynam- 
ical system representing the entire set of economic inter- 
actions in a society, analogous to those previously derived 
for gases and liquids. By viewing the economy as a ther- 
modynamic system, one can identify the income distribu- 
tion with the distribution of energy among the particles 
in a gas. In particular, a class of kinetic exchange models 
have provided a simple mechanism for understanding the 
unequal accumulation of assets. Many of these models, 
while simple from the perspective of economics, has the 
benefit of coming to grips with the key factor in socioeco- 
nomic interactions that results in very different societies 
converging to similar forms of unequal distribution of re- 
sources (see Refs. 0,i], which consists of a collection of 
large number of technical papers in this field). 

In this paper, we consider a new model of asset ex- 
change on a directed network and investigate the nature 
of wealth distribution, using numerical simulations. In 
Sec. [II] we review the salient features of the gas-like ki- 
netic models. In Sec. IIII1 we propose the new model 
and present the results. We conclude with discussions in 
Sec. HV] 



II. GAS-LIKE MODELS 

In analogy to two-particle collision process which re- 
sults in a change in their individual kinetic energy or 
momenta, income exchange models may be defined using 
two- agent interactions: two randomly selected agents ex- 
change money by some pre-defined mechanism. Assum- 
ing the exchange process does not depend on previous 
exchanges, the dynamics follows a Markovian process: 



mj{t 



= M 



m,i(t) 
mj(t) 



(2) 
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where mj(i) is the income of (individual or corporate) 
agent i at time t and the collision matrix M defines the 
exchange mechanism. 

In this class of models, one considers a closed economic 
system where the total money M and number of agents N 



2 



is fixed. This corresponds to no production or migration 
in the system where the only economic activity is con- 
fined to trading. In an y tr ading, a pair of traders i and j 
exchange their money jl6l . UtI Ha . [l9 ] , such that their to- 
tal money is locally conserved and nobody ends up with 
negative money (rrii(t) > 0, i.e, debt not allowed): 

rrii(t + 1) = nii{t) + Am; rrij(t + 1) = Tnj(t) — Am. (3) 

Time (t) changes by one unit after each trading. 

The simplest model considers a random fraction of 
total money to be shared [13]. At steady-state (t — > 
oo) money follows a Gibbs distribution: P(m) = 
(l/T)cxp(-m/T); T = M/N, independent of the ini- 
tial distribution. This follows from the conservation of 
money and additivity of entropy: 



XjUijit)- 



P(mi)P(ni2) = P(mi + m 2 ). 



(4) 



This result is quite robust and is independent of the 
topology of the (undirected) exchange space, be it regular 
lattice, fractal or small- world [2(| ■ 

A saving propensity factor A was introduced in the 
random exchange model [l9| , where each trader at time t 
saves a fraction A of its money rrii(t) and trades randomly 
with the rest: 

rrii(t + 1) = Xm^t) + e y - [(1 - A)(m i (t) + m,j(t))] , (5) 



m 3 (t + 1) = Xm 3 (t) + (1 - ey) [(1 - \)(rm(t) + mj(t))] , 

(6) 

eij being a random fraction, coming from the stochastic 
nature of the trading. 

In this model (CC model hereafter), the steady state 
distribution P(jn) of money is decaying on both sides 
with the most-probable money per agent shifting away 
from m = (for A = 0) to M/N as A -> 1 [11. This 
model has been understood to a certain exte nt 1221 l23j 
and argued to resemble a Gamma distribution 22]. But, 
the actual form of the distribution for this model still 
remains to be found out. It seems that a very similar 
model was proposed by Angle [U [25| several years back 
in sociology journals. The numerical simulation results 
of Angle's model fit well to Gamma distributions. 

Empirical observations in homogeneous groups of in- 
dividuals as in waged income of factory labourers in UK 
and USA Q and data from population survey in USA 
among students of different school and colleges produce 
similar distributions (25j |. This is a simple case where a 
homogeneous population (say, characterized by a unique 
value of A) has been identified. 

In a real society or economy, saving A is a very inho- 
mogeneous parameter. The evolution of money in a cor- 
responding trading model (CCM model hereafter) can be 
written as (26j |: 

mj(t+l) = X l m l {t) + e lj [(1 - X l )m i {t) + (1 - Xj)mj(t)] , 



-e y -)[(l-A l KW + (l-A,K-W] 

(8) 

The trading rules are same as CC model, except that Aj 
and Xj, the saving propensities of agents i and j, are dif- 
ferent. The agents have fixed (over time) saving propen- 
sities, distributed independently, randomly as A(A), such 
that A(A) is non- vanishing as A — > 1, Ai value is quenched 
for each agent (Ai are independent of trading or t). The 
actual asset distribution P(m) in such a model will de- 
pend on the form of A(A), but for all of them the asymp- 
totic form of the distribution will become Pareto-like: 
P(m) ~ m~( 1+v >; v = 1 for m — » oo. This is valid 
for all such distributions, unless A(A) oc (1 — A)" 5 , when 
P(m) ~ m-( 2+s ~> [II, [III lH. In the CCM model, agents 
with higher saving propensity tend to hold higher aver- 
age wealth, which is justified by the fact that the saving 
propensity in the rich population is always high (29| . An- 
alytical understanding of CCM model has been possible 
until now under certain approximations [30l ]. and mean- 
field theory [27j • Unlike the above models with savings 
as a quenched disorder, one can also consider savings as 
an annealed variable and still derive a power law distri- 
bution in wealth 1311. 



III. MODELS ON DIRECTED NETWORKS 

The topology of exchange space in a real society is 
quite complicated. A mean field scenario does not take 
into account the constraints on the flow of money or 
wealth. In other words, there are strong notions of di- 
rectionality and sometimes, hierarchy in the underlying 
network, where money is preferentially transeferred in 
certain direction that others, contributing in irreversible 
flow of money. A way to imitate this is to consider wealth 
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FIG. 1: The distribution R(pi) of the link disorder 
Pi for a directed network with different values of p = 
0.50, 0.30, 0.10, 0.05, 0.01 for a system of TV = 100 nodes, ob- 
tained by numerical simulation, averaged over 10 4 realiza- 
tions. 
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exchange models on a directed network [H, [H, l34j . 
There have been previous attempts to obtain the same 
using the physics of networks [35LI36L H3] ■ In this section, 
we would propose suitable toy models and present some 
interesting observations. 

A. The Model 

We consider N agents, each at a separate node, each 
of which are connected to the rest N — 1 by directed 
links. The directionality of the links denote the direction 
of flow of wealth in this fully connected network. The 
directed network parametrized by p is constructed in the 
following way: 

(a) there are no self links, so that the adjacency matrix 
A HUHl] has diagonal elements a,, = for all sites i. 

(b) for each matrix element , i ^ j, we call a random 
number r. ; = +1 if r < p and = — 1 otherwise. 
Also a,ji = —dij. Thus we have N(N — l)/2 such calls 
of r. ciij — +1 denotes a link directed from i to j and 
fly = — 1 denotes a link directed from j to i. 

We define pi = j^zi a u as a measure of link dis- 
order at the site i. • denotes the sum over all N — 1 
sites j linked to i. Thus, pi — 1 is a node for which all 
links are outgoing and pi = — 1 is a node for which all 
links are incoming. The parameter p has a symmetry 
about 0.5 and the distribution R(pi) is also symmetric 
about 0, which is, in fact, a consequence of the conser- 
vation of the number of incoming and outgoing links. A 
network with p = 0.5 has the lowest degree of disorder, 
given by a narrow distribution R{p) of p, around p = 
(see Fig.[IJ. This means that almost all nodes have more 
or less equal number of incoming and outgoing links. On 
the other extreme, p = 0.01 is a network which has a 
small but finite number of nodes where most links are 
incoming/outgoing, thus giving rise to a very wide dis- 
tribution of link disorder R(p). The rules of exchange are 
as follows: If ay = — 1, 

rriiit + 1) = rriiit) + pjrrij{t) 
mj(t+l) = mj(t) - Hjirijit) 

(9) 

else, ciij = +1, 

rm(t + l) = mi(t) — mmi{t) 
rrij(t + l) = mj(t) + pinii(t). 

(10) 

where < pi < 1 is a 'transfer fraction' associated with 
the ith agent, and rrii(t) is the money of agent i (or, 
money at node i) at time t. The total money in the 
system is conserved, no money is created or destroyed, 
as is evident from eqn. (J9]) and eqn. (fT0|) . If there is a 
link from j to i, the node i gains pj fraction of jth agent's 
money, which, of course agent j loses. Otherwise, if there 
is a link from i to j, the node j gains pi fraction of ith 



agent's money, which, of course agent i loses. In the 
Monte Carlo simulations, one assigns random amount of 
money to agents to start with, such that the average 
money M/N =1. A pair of agents (nodes) are chosen 
at random, and depending on the directionality of the 
link between them (the sign of dij), the relevant rule, 
eqn. © or eqn. (jTUJ) is chosen. This is repeated until 
a steady state is reached and the distribution of money 
does not change in time. The distribution of money P(m) 
is obtained by averaging over several ensembles (different 
random initial distribution of money). 

This model is different from the CC and CCM models, 
but one can relate the transfer fraction p analogous to A 
in the CC and CCM models. 



B. Model with uniform p 

First we discuss the case of homogeneous agents, i.e, 
when all agents i have pi = p. The p = limit is trivial, 
as the system does not perform any dynamics. Fig. [^a) 
shows the steady state distribution P{m) of money m 
for p — 0.1 for different values of network disorder p. 
In general, the distribution of money has a most proba- 
ble value, which shifts monotonically from about 0.85 for 
p = 0.5 to as p — > 0. P{m) has an exponential tail, 
but a power law region develops as p — > in between the 
most probable peak and the exponential cut-off, which 
fits approximately to m~ 15 . For the limit p — ► 0, con- 
densation of wealth at the node(s) with strong disorder 
(p — > 1) is apparent from the single, isolated data point at 
the m max = M end (see inset of Fig. [Ha), for p = 0.01). 
There is a strong finite size effect involved in this be- 
havior. To emphasize this, we plot P{m) for p = 0.01 
for N = 100,500 and 1000 (inset of Fig. [S^a)). While 
N = 100 and N = 500 does show the isolated data point, 
it is absent for N = 1000. This indicates that this behav- 
ior is absent for infinite systems for p — > 0. For larger val- 
ues of p, the distribution resembles Gamma distributions, 
as in the CC model. For p = 0.5, the most-probable value 
of P(m) is always at (see Fig. [2Kb)). For weak disorder 
(p = 0.5), P(m) is exponential, but it goes to a wider 
distribution as one goes to higher disorder (p — > 0). The 
condensation of wealth at node(s) with high value of p 
(p — > 1) is again apparent from the single, isolated data 
point at the m max = M end (see Fig. HJb), for p = 0.05). 
For p — 0.9, P(m) is always decaying, with a wide distri- 
bution upto m max = M (see Fig.[^c)). As like previous 
plots, the condensation of wealth at node(s) with high 
value of p is visible: see plot for p = 0.05 in Fig. (D(c). A 
common feature for the curves for all values of p is that, 
P(m) exhibits log-periodic oscillations, while resembling 
roughly a power-law decay. Another important feature is 
that P(m) — > N for m — > 0, which indicates that money 
is distributed in a very small fraction of nodes, while most 
nodes have almost no money at a given instance. 

For a particular value of the network disorder p, the 
wealth distribution P(m) becomes more and more 'fat 
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FIG. 2: The steady state distribution P(m) of money m for directed networks characterized by different values of p. (a) For 
/i = 0.1, the inset shows P(m) for p = 0.01 for N = 100, 500, 1000, and the power law m~ 1,5 is also indicated, (b) For /i = 0.5 
and (c) /i = 0.9. (d) shows the plots for uniform, random distributed /i, D(fi) — 1, and also a guide to the power law m~ 2 . 
The data are obtained by numerical simulation, for a system of N = 100 nodes, averaged over 10 3 realizations in the steady 
state and over 10 4 initial configurations. The average money M/N is 1. 



tailed' as /i is increased. This is in contrast to what is 
observed in CC model where P(m) organizes to a 
narrower distribution as A increases. 



C. Model with distributed /i 

We now consider the case when agents i have have dif- 
ferent values of /Ltj, which do not change in time. This 
is a case of heterogeneous agents where the heterogene- 
ity can be viewed as a 'quenched disorder'. We consider 
a random uniform distribution of /i, i.e, D(jj,) = 1 in 
< H < 1. This is the case analogous to the CCM 
model [2q . Fig. [21(d) shows the plots of the money dis- 
tribution P(m) for different values of network disorder p. 
All curves have a power law tail, resembling a m~ 2 vari- 
ation. However, the effect of the topology of the under- 
lying network are visible: For strong disorder in topol- 
ogy P — 0.05, condensation of wealth at node(s) with 
high value of p (p — ► 1) is also apparent from the sin- 



gle, isolated data point at the vn max = M end. Further 
investigations also indicate that the power law exponent 
is similarly related to the distribution of 'transfer frac- 
tion' //, as one observes in the CCM model [H, i-e, 
P(m) ~ m" 2 for most distributions, while one can obtain 
P[m) ~ m-( 2+l5 ) if D{p) oc (1 - ^i) s . 



IV. DISCUSSIONS 

The empirical data for income and wealth distribution 
in many countries are now available, and they reflect a 
particular robust pattern: The bulk (about 90%) of the 
distribution resemble the century-old Gibbs distribution 
of energy for an ideal gas, while there are evidences of 
considerable deviation in the low income as well as high 
income ranges. The high income range data (for 5-10% 
of the population in any country) fits to a power law tail, 
known after Pareto, and the value of the (power law) ex- 
ponent ranges between 1-3 and depends on the individual 
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make-up of the economy of the society or country. 

The analogy with a gas like many-body system has 
led to the formulation of the models of markets. The 
random scattering-like dynamics of money (and wealth) 
in a closed trading market, in analogy with energy 
conserved exchange models, reveals interesting features. 
Self-organization is a key emerging feature of these ki- 
netic exchange models when saving factors are intro- 
duced. These models [ljl [26[ produce asset distribu- 
tions resembling that observed in reality, and are quite 
well studied now [28| . Empirical observations in homoge- 
neous groups of individuals as in waged income of factory 
labourers in UK and USA [7|] and data from population 
survey in USA among students of different school and col- 
leges produce similar distributions [25[. There have also 
been other models [38] which uses the kinetic exchange 
approach. 

Heterogeneity is an inherent feature of real networked 
economy, mostly observed with with behavior of agents 
contributing to irreversible flow of money. A simple way 
of looking into it is to consider heterogeneity in the in- 
teraction topology for the agents. In all earlier studies, 
the effect of the topology of the space on which the ki- 
netic exchanges took place, was overlooked. Most of the 
models were either on regular, undirected lattices, or on 
undirected graphs. Our study concentrates on the effect 
of a directed network as the exchange space for such ki- 
netic exchanges. We consider a simple, fully connected, 
directed network. A notion of disorder p in the directed- 
ness of the graph was introduced, in analogy to real world 
networks, which depend on the number of incoming and 
outgoing links at a particular node. The distribution 
R(p) characterizes a particular network. In analogy to 
saving propensity A in CC and CCM models [H, [26| , we 
define a 'transfer fraction' /i. We consider two distinct 
models: (i) with uniform p, (Sec. IIII B[) . and (ii) with 
distributed p fSec. IIII"C|) . We employ numerical simula- 
tions to study the above models. For uniform p, when the 
transfer fraction p is small, low degree of disorder in the 
network produced money distribution P(m) with most- 
probable money at finite values of m similar to Gamma 
distribution observed in CC model [lj|, while high dis- 
order produced a wider distribution of P(m) (Fig.[2{a)). 
At intermediate ranges of p, P(m) has its most-probable 
value always at m — > (Fig. H[b)), while for high values 
of (i, P(m) has a wide distribution with log-periodic os- 
cillations (Fig. HIc)). For the model with p distributed 
randomly and uniformly (D(p) — 1), money distribu- 
tion exhibits power law P(m) ~ m~ 2 for large money m 
(Fig. [U(d)). The common feature of all these cases was a 
very wide distribution P(m) of money, when the network 
disorder is very strong (see plots for p = 0.05 in Fig. [2]), 
accompanied by a signature of accumulation of money 
at the nodes/with the agents having strongest disorder, 



which is observed in the data point at m m ax — M. How- 
ever, this effect is more pronounced in smaller system 
sizes, while it may disappear in the thermodynamic limit 
(see inset of Fig. [2£a)). 

A simple model of similar exchanges have been studied 
before [T3|. The case corresponds to the choice p = 0.5 
for our model. Moreover, the transferred money was in- 
dependent of m and not proportional to p. Money distri- 
bution was found to be exponential despite explicit viola- 
tion of time-reversal symmetry. In an earlier paper [13], 
results for multiplicative random exchanges correspond 
to certain cases of the present paper. It was found that 
P(m) vanishes at m — for p < 0.5 and diverges for 
p > 0.5, and is exponential for p = 0.5. These results 
correspond to that of our model for p = 0.5 (Sec. IIIIBj) . 

We have mostly used the terms 'money' and 'wealth' 
interchangeably, treating the models in terms of only one 
quantity, namely 'money' that is exchanged. Of course, 
wealth does not comprise of (paper) money onlv, and 
there have been studies distinguishing these two [28|, [3^, 

Study of such simple models here give some insight 
into the possible emergence of self organizations in such 
markets, evolution of the steady state distribution, emer- 
gence of Gamma-like distribution for the bulk and of the 
power law tail, as in the empirically observed distribu- 
tions. The study of models on directed graphs gives us 
new insights: (i) one can obtain wide asset distribution 
even in homogeneous populations (as in uniform p), (ii) 
the power law exponent for the distributed p case behaves 
in the same manner as in CCM model, and (iii) strong 
disorder in topology of the underlying directed network 
produces accumulation (at p — > 1) or the opposite effect 
(at p — ► — 1) of wealth, thus giving rise to a higher degree 
of wealth inequality in the system. 

These model studies also indicate the appearance of 
self-organization, and the self-organized criticality [4l| in 
particular, in the simple models. These have prospective 
applications in other spheres of social science, as in appli- 
cation in policy making and taxation, and also physical 
sciences, in designing desired energy spectrum for differ- 
ent types of chemical reactions [42[. 
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